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Abstract 

We investigate chemical-potential (fi) and temperature (T) dependence of scalar and pseudo-scalar me- 
son masses for both real and imaginary /x, using the Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) 
model with three-flavor quarks. A three-flavor phase diagram is drawn in fi'^-T plane where positive (neg- 
ative) /i^ corresponds to positive (imaginary) /x. A critical surface is plotted as a function of light- and 
strange-quark current mass and fi^. We show that /i-dependence of the six-quark Kobayashi-Maskawa-'t 
Hooft (KMT) determinant interaction originated in C/a(1) anomaly can be determined from lattice QCD 
data on t]' meson mass around /x = and fi = ittT/S with T slightly above the critical temperature at 
/i = where the chiral symmetry is restored at /x = but broken at /x = ivrT/S, if it is measured in future. 
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I. INTRODUCTION 



In recent theoretical studies, novel scenarios for QCD phase structure at finite real chemical 
potential (/xr) are suggested; for example, the quarkyonic phase i^iM, the multi critical-endpoint 
generation [15|48|] and the Lifshitz-point induced by the inhomogeneous phase [|9|]. Thus, qualitative 
or speculative investigation of QCD phase diagram is progressing well. 

Nevertheless, quantitative or more conclusive understanding of QCD phase diagram is quite 
poor. The principal reason is the sign problem in the first-principle lattice QCD simulation at finite 
yUR. Several methods such as the reweighting method fllQI . the Taylor expansion method [lllll . the 
analytic continuation from imaginary chemical potential fii to /iR [Il2l-ll4ll and so on were proposed 
so far to circumvent the sign problem. However, they do not reach the /Ur/T > 1 region yet. For 
this reason, effective models such as the Nambu-Jona-Lasinio (NIL) model were used so far to 
investigate qualitative properties of the phase structure at finite /xr. The effective-model approach, 
however, has an ambiguity particularly in the interaction part; see Ref. [llSll and references therein. 

Thus, a new approach should be proposed for quantitative or more reliable investigation of QCD 
phase diagram at finite /iR. As a possible answer, recent ly, vye p roposed the imaginary chemical 



potential matching approach (the /ij-matching approach) [|16l ll7ll. In this approach, interactions of 



the effective model are determined from LQCD data at finite jii where no sign problem comes out. 
After the determination, a phase structure at finite /xr is predicted with the effective model. The 
most important point in this approach is whether the model taken can reproduce the Roberge- Weiss 
(RW) periodicity and the RW transition at finite /xi [18]. In our previous works ll^, we showed 



that the Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) model [|l9ll can do it, because the 
thermodynamical potential of the PNJL model is invariant under the extended Z3 transformation 



of 



-)■ e 



±ie+i- 



where 9 = Hi/T. Here, <1> and ^ denote the Polyakov-loop and its conjugate, respectively. This 
symmetry ensures the RW periodicity. Since, the PNJL model is designed to treat the confinement 
mechanism approximately in addition to the chiral symmetry breaking, we can investigate not 
only the chiral transition but also the deconfinement transition with the PNJL model. We also 
showed by using the PNJL model that the crossover deconfinement transition that takes place at 
finite 9 becomes stronger as 9 increases and eventually at 6* = 7r/3 it changes into the RW phase 
transition nluA . 



; 2-Kk 
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The Ua{1) anomaly related to instantons can be taken into account in the NJL and PNJL mod- 
els. In the three-flavor case, it is described by the effective six-quark Kobayashi-Maskawa-'t Hooft 



(KMT) determinant interaction [|2(a 



2l\\ . The t/A(l) anomaly restoration at finite T in the case of 



fi = is investigated by the NJL model [|22l1 that can reproduce the lattice QCD data [|23|]. At 
finite /iR, /i-dependence of the anomaly restoration strongly depends on that of coupling constant 
of the KMT determinant interaction \\24i\ . However, yU-dependence of Gd is unclear, because 
lattice QCD data is not feasible at finite /xr and also theoretical understanding on yU-dependence 
of the instanton density is not sufficient. Therefore, the phase structure in the three-flavor system 
is more ambiguous than in the two-flavor system. 

In this paper, we investigate scalar and pseudo-scalar meson masses in both the /ir and /ij 
regions, using the three-flavor PNJL model. We show i]' meson mass is sensitive to Go particularly 
near 9 = n/3. This means that the 9 dependence of r/' meson mass is a good quantity to determine 
fi dependence of Go- At the present stage, there is no reliable lattice QCD data particular on 
meson masses for the case of finite j^i. Therefore, our investigation is limited to only a qualitative 
level. 



II. THREE-FLAVOR PNJL MODEL 

Lagrangian density of the three-flavor PNJL model is 

8 

-CpNjL =q{iiuD'' - mo)q + Gs "^[{qKqy + {qiibKq) 



[\q^aq) ^ yqfjb^aq)'^ 

a=0 

— Gd 



detg,(l + 75)g, + detg,(l-75)g, - T), (2) 

L ij tj 



where D'^ = d" + iA" = d'' + iS^gA^^Xa/^ with the gauge coupling g and the Gell-Mann matrices 
\a- Three-flavor quark fields q = {qu, qd, qs) have current quark masses mo = diag(m„, m^, m^). 
The Polyakov potential U is defined later in (|7]) and ([8]). In the interaction part, Gs and Go denote 
coupling constants of the scalar-type four-quark and the KMT determinant interaction, respec- 
tively. The determinant det runs in the flavor space and then the KMT determinant interaction 

ij 

breaks the Ua{1) symmetry explicitly. 

In the PNJL model, the gauge field is treated as a homogeneous and static background field. 
The Polyakov-loop ^ and its conjugate ^ are given by 

^ = ^tr,(L), <P = itre(L) (3) 
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where L = exp{iA4/T) with A4 = iAo in Euclidean space. In the Polyakov-gauge, is diagonal 
in the color space. 

We make the mean field approximation (MFA) to the quark-quark interactions in Q in the 
following way. In the operator product q^qj is first divided into q^qj = aij + (qiqj)' with the 
mean field cTj^ = (qiqj) and the fluctuation (qiqj)' where i, j = u,d,s. Ignoring higher-order terms 
of (qiqj)' in the rewritten Lagrangian and re-substituting {qiqj)' = qiqj — (Tij into the approximated 
Lagrangian, one can obtain a linearized Lagrangian based on MFA: 

i=u,d,s 

-W(<?[A],^[A],T), (4) 

where the dynamical quark mass Ma is defined by Ma = rrii — AGscth + 2G£)0'jjO'kk with % 3 
k. 

In this study, we impose the isospin symmetry for u-d sector and then we use mi = rriu = md. 
The thermodynamical potential becomes 



f — J „ ^ 



f=u,d,: 



(27r) 



NrE„ 



+ ( ^Gsal - AGoauu^^ddCTss) + U{<P[A],^A],T). (5) 



i=u,d,i 

We take the three-dimensional momentum cutoff, 



"^'^ tdpp\ (6) 



(27r)3 27r 

because this model is non-renormalizable. Hence, the present model has five parameters Gs, Gd, 
mi, ms and A. We use the parameter set determined in Ref. jlS]; these are fitted to the empirical 
values of vr meson mass and its decay constant, K meson mass and its decay constant and y' meson 



mass. We also use U of Ref. [,26,1 fitted to LQCD data in the pure gauge limit at finite T [|27l . 



281]: 



h{T) = ao + a, (^) + a, (^) ' + as (^)'. (8) 



In this study, we take the original value To = 270 MeV. 



III. MESON MASS FORMALISM 



First, Lagrangian density Q is rewritten by MFA into 



a=0 



+ X] abi(l^aq){qKq) + Gl^^{qi-i^\aq){qil^\bq)\ 

a,b=0,3,8 



(9) 



with 



Gq — Gs =F -;;Gd(o'uu + <^dd + 



— G2 — ^3 — Gs ± -Gof^s 



4 — ^5 
1 



Gs ± ■^^GoO'dd, 



G^ 



Gf — Gs ± -GD<yu 



Cg — Gs ± -Goi'^Ouu + 2arfrf — cr^^), G3Q — Gq3 — =F — 1=Gd{Ouu — crdd), 




G, 



08 ~ '-^80 



V2 ± 

±Y^Gz)(cr„„ + (Jdd — 2(Tss), Ggg = G, 



± 

83 



1 

1 

"273 



(10) 



GnicTuu — CTdd)- 



In the present case that the isospin symmetry is imposed in the u-d sector, we have G^g = G 



03 



^38 — ^83 — 0- 



Taking the same procedure as in the two-flavor case [|l7h . we can obtain dynamical meson 
masses from poles of the effective propagator ll25ll : 

22G,- 



^ 1 - 2G,i7(go; 







qo=M(_ 



(11) 



1 - 2Gii7(go) 

where U denotes the Polarization function. The subscript i stands for meson ,^ in a state i. The 
polarization function between states i and j is represented by 



rTr 



r{T,)SF{p)r{Tj)Sp{p-q) 



(12) 



where Sp denotes the quark propagator and the vertex function r is 1 for the scalar meson and 
^75 for the pseudo-scalar meson. The matrices Tj in flavor space depend on meson considered; 
for example, Tj = Tj = A3 for vr and mesons and Tj = (Ae + iX7)/V^,Tj = (Ae + iX-j)!^ 
for K and k mesons. When T and are finite, the corresponding equations are obtained by the 
replacement 



Po i^n + fJ' ~ = iT^T{2n + 1) + yU — Zy44, 



(2^ 



d^p 
(2^ 



(13) 
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As for Tj and rj' mesons, the effective coupling constant and the polarization function U are 
2x2 matrices in flavor space, since the isospin symmetry is imposed for the u-d sector: 



'-^80 '-^88 



n 



(14) 



where 



iloo = i [277" + 77^^] , TTss = i [TT" + 27r^^] , TTos = TTgo = ^ [TT" - TT^^] . (15) 



Here, the polarization function TT^-^ for each flavor / is defined by 



Tr 



C,D 



rsUp)rsUp-q) 



(16) 



with Sp{q) the propagator of quark with flavor /, where in Ttc,d the trace is taken for color and 
Dirac indices. Therefore, r] and t]' meson masses satisfy 



det(l - 2G+7Z") = 0. 



(17) 



This equation has two solutions; the lower corresponds to r] meson mass, while the higher does to 
1]' meson mass. Masses of a and /o meson are obtained by replacing (7+ by G^ in (flTl) and setting 
r = 1 in (fill). 



IV. NUMERICAL RESULTS 

Figure [T] shows the phase diagram in /i^ — T plane for the three-flavor case. The positive 
(negative) fj? half-plane means real (imaginary) /i. In the present parameter set, a critical endpoint 
arises in the positive /i^ half -plane, while a RW endpoint does at /i^ = — (ttT/S)^ with T ^ 
250 MeV in the negative /i^ half -plane. At the RW endpoint, the phase transition is second order. 



In the two-flavor case, as shown in [291], the order of RW phase transition at RW endpoint depends 



on the Polyakov potential U taken; it is second order for U of Ref. [lol], but first order for U of 
Ref. [30]. And the latter gives a result more consistent with lattice QCD data at finite jii than the 
former. This sort of analysis is quite important also for the three-flavor case in order to determine 
the form of U, if precise lattice QCD data on the RW endpoint become available in future. 

The critical endpoint (closed circle) in Fig. [U is a function of mi, and /x^. This is described 
as a surface in the nii-nis-^'^ space. The surface, usually called the critical surface, is plotted 
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Fig. 1: (Color online) Phase diagram in ^u^ — T plane for the three-flavor case. The solid (dotted) line 
denotes a first-order (crossover) chiral transition, while the dot-dashed (dashed) line does a first-order RW 
(crossover deconfinement) transition. The closed circle (square) stands for an endpoint of the first-order 
chiral (RW) phase transition. 
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Fig. 2: (Color online) Critical surface as a function of mi, rus and //^ in the //r and /ij regions. 

in Fig. [21 In the present model, the critical surface has a positive curvature and then a critical 
endpoint arises somewhere in the fi^-T plane, as shown in Fig. [U when mi and are taken to 
be physical values. This result may change, if the coupling constant Gd of KMT determinant 
interaction depends on /x [.24.1 . However, the /x-dependent Gz) of Ref . jl^ can not be applicable to 
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the imaginary /i region, since it breaks the extended Z3 symmetry, i.e. the RW periodicity. Thus, 
at the present stage, we have no way of determining jj, dependence oiGo- We then do not consider 
any ^u-dependent G/) in this paper. 

In Fig.[3l we investigate yU^-dependence of meson mass for mesons (tt, K, r], r]', a, k, oq, /o), 
using the PNJL model in which Go is a constant Go = G^iO) = —12.36/1^^ where G'd(O) 
is determined at T = /i = 0. The left and the right panels correspond to T = 200 and 300 
MeV, respectively. Furthermore, Figure |4] shows 9 dependence of meson mass for T = 300 MeV; 
the left (right) panel corresponds to tt and a (K and k). Obviously, these meson masses are 9- 
even and have the RW periodicity. Pion mass has dips at ji^ ^ 0.06 (GeV)^ in panel (b), at 
/i^ ^ -0.08 (GeV)2 in panel (d) of Fig.[3]and at ^ ^ ±0.7 and ibl.3 MeV of Fig. H These are 
threshold effects due to tt — )■ quark + antiquark. 

A mass difference between the chiral partners ir and cr is a good indicator of the chiral symme- 
try; the symmetry is restored (broken) when the difference is small (large). As for T = 300 MeV, 
as shown in panel (c) of Fig. [3] and the left panel of Fig. |4l the chiral symmetry is restored at 
/i^ > -0.02 (GeV)2, but broken at fi^ < -0.02 (GeV)^ (vr/G <9< n/3). As for T = 200 MeV, 
as shown in panel (a) of Fig. [3l the symmetry is restored at fi^ > 0.1 (GeV)^, but broken at 
yU^ < 0.1 (GeV)^. As an interesting result in Fig.[3l panels (c) and (d) almost agree with panels (a) 
and (b), respectively, if in panels (c) and (d) the fi'^ scale is shifted to the left by about 0.1. Thus, 
shifting the yU^ scale to the left corresponds to looking at meson mass at lower T. If meson mass 
is measured in future by lattice QCD in the negative //^ region for some temperature TLatt^ the 
behavior qualitatively agrees with that in the positive /x^ region for temperature lower than TLatt- 
Therefore, we can predict qualitative behavior of meson mass in the positive fi^ region from lattice 
data on meson mass in the negative region. 

Next, it is investigated how the KMT determinant interaction affects meson mass at imaginary fi 
by changing the value of Gd from the original one GoiO)- As mentioned above, theoretically, Gd 
is allowed to have fi dependence. However, since the actual form is unknown, we simply change 
the value of Gd in the present analysis. Figure [5] shows tt and a meson masses at T = 300 MeV 
with the KMT determinant interaction in which Gd = for left panel and Gd = G'£)(0)/2 for 
right panel. AsGd increases, the tt meson mass is reduced with almost keeping the 9 dependence. 
The cr meson mass is also reduced, but the 9 dependence is changed a lot around 9 = ir/S. The 
left (right) panel of Fig. [6] shows 1] (r/') meson mass at T = 300 MeV in three cases of Gd = 0, 
Gd{^)/'^ and Gd{^)- The rj meson mass has a similar property to the cr meson mass. Most 



8 




Fig. 3: (Color online) /x^-dependence of meson mass for (a) tt, a, K and k at T = 200 MeV, (b) 1], 1]' , oq 
and /o at T = 200 MeV, (c) vr, a, K and k at T = 300 MeV, and (d) r?, r/', oo and /o at T = 300 MeV. 

interesting and important property is 6 dependence of rj' meson mass. The rf' meson mass has a 
weak Gd dependence at = 0, but the Gd dependence becomes strong around = n/3. Thus, 9 
dependence of Gd that is allowed theoretically can be determined from that of r]' meson mass, if 
it is measured in future by lattice QCD. 

The sensitivity of i]' meson mass to Go around 9 = n/3 can be understood in the following. 
The KMT determinant interaction affects the meson mass only through the dynamical quark mass 
Mjj = rrii — AG s(Tii + 2G £,a jjakk with a G/j dependent term of form GoO'iiO'jj- The term GDCuO'jj 
is strongly suppressed when the chiral symmetry is restored, even if is large. Hence, the f/A(l) 
anomaly affects the meson mass only when the chiral symmetry is broken. As shown in Fig.[Tl the 
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Fig. 5: (Color online) ^-dependence of a and tt meson masses at T = 300 MeV for Gd = in the left 
panel and for G/j = G£i(0)/2 in the right panel. 

critical temperature of the chiral transition goes up as 6 increases from to vr/S. Thus, 77' meson 
mass is most sensitive toG'£)at^^ = 7r/3 when 6 increases from to vr/S with temperature fixed, 
because an and ajj are largest there. 

In this study, is assumed to be constant. If rj' meson mass is measured by three-flavor lattice 
QCD in future, there is a possibility that the PNJL model with constant Gd can not reproduce lat- 
tice QCD data. If so, the deviation can determine 6'-dependence of and hence yUR-dependence 
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Fig. 6: (Color online) 6'-dependence of meson mass at T = 300 MeV for rj in left panel and -q' in right 
panel. Three values of = 0, Gi:)(0)/2 and Gd{0) are taken. 

of Gd. 

V. SUMMARY 

Using the three-flavor PNJL model, we have analyzed yU-dependence of scalar and pseudo- 
scalar meson masses in both the real and the imaginary /i region. In the imaginary fi region, 
the meson masses are even functions of 9 with the RW periodicity. As an interesting result, /i^ 
dependence of meson mass in the negative region at some temperature is close to that in the 
positive /i^ region at temperatures lower than the temperature. If meson mass is measured in future 
by lattice QCD in the negative fi^ region for some temperature TLatt, the behavior qualitatively 
agrees with that in the positive fi'^ region for temperatures lower than TLatt- Therefore, we can 
predict qualitative behavior of meson mass in the positive fi'^ region from lattice data on meson 
mass in the negative /i^ region. 

The t/A(l) anomaly (the KMT determinant interaction) affects meson masses through the term 
GoCTjjakk in the dynamical quark mass Ma. Particularly, the effect is remarkable for r]' meson. 
For temperatures slightly above the critical temperature at = 0, the chiral condensate increases a 
lot as 9 increases from to tt/S, so that the effect has a strong 9 dependence. We then recommend 
that meson masses, particularly 7]' meson mass, be measured by lattice QCD for = and n/3 at 
such higher temperatures. Using the lattice QCD data, we can determine T and fi dependences of 
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coupling constant Go of the KMT determinant interaction and hence can predict the three-flavor 
phase diagram with higher reliability by the PNJL model. 
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